#include <buola/mat.h>
#include <buola/mat/prettyprint.h>
#include <buola/mat/svd.h>
#include <buola/mat/traits_typename.h>

using namespace buola;

int main(int /*argc*/,char** argv)
{
    mat::CMat_d a=mat::random(2000,2000,std::normal_distribution<double>());
    mat::CMat_d b=mat::random(2000,2000,std::normal_distribution<double>());

    start_timer();
    double lSum;
    for(int i=0;i<100;i++)
    {
        mat::CMat_d c=(a+2*b)/(3*(b+a));
        lSum+=c(3,3);
    }
    end_timer();
    msg_info() << lSum << "\n";
}